import sympy as sym
def eq(x,y):
    return sym.Eq(x,y)

f = sym.symbols("f")
print(sym.solveset(sym.Eq(f**2,1),f))
x = sym.symbols("x")
tan = sym.sin(x) / sym.cos(x)
q = eq(tan,sym.tan(x))
print(q)
